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ABSTRACT 

The fastest millisecond pulsar PSR B1937+21 presents an interpulse separated from the main pulse 
by nealy 180° at radio frequencies. Recently, the ASCA observations (Takahashi et al. 2001) detected 
pulsed X-ray emission from this pulsar. Only a single narrow X-ray pulse is observed, which is coincident 
with the radio interpulse in phase. We investigate the possible origin of the pulsed X-rays from the polar 
cap (PC) accelerators or the outer gap (OG) accelerators in the frame of PC model and OG model, 
respectively, by assuming a dipolar magnetic field structure and the same radio emission pattern from 
its poles for the pulsar. The OG model can naturally explain the main observational facts. For the PC 
model, the coincidence between the X-ray pulse and the radio interpulse can not be reproduced in the 
assumed case. However when considering possible deviation from our assumption, PC model may still 
be valid for this pulsar in some cases. 

Subject headings: pulsar: PSR B1937-t-21: X-rays: emission mechanism 



1. INTRODUCTION 

Even more than 30 years after the discovering of high 
energy pulsars, the theoretical reproduction of X-ray and 
7-ray emission from such pulsars is still a matter of de- 
bate. It is commonly agreed that there are two scenar- 
ios on modelling X-ray and 7-ray creation: the outer gap 
(OG) model (e.g., Cheng, Ho & Ruderman 1986a, b, Ro- 
mani 1996, Cheng & Zhang 1999, Hirotani & Shibata 
2001) and the polar cap (PC) model (e.g., Harding 1981, 
Sturner & Dermer 1994, Daugherty & Harding 1996, Luo 
et al. 2000, Zhang & Harding 2000). The fundamental dif- 
ference between these two types of models is the location 
of the regions where particles are accelerated to relativis- 
tic energies and emit high energy photons. The early PC 
model (Harding 1981) assumed that the emission is pro- 
duced just above the PC surface. In the present versions 
of PC model (e.g. Daugherty & Harding 1996, Harding 
& Muslimov 1998), it is proposed that the particle accel- 
eration region may extend from the PC surface to several 
stellar radius above due to free charge flow and inertial 
frame dragging, so that wide double-peak X-ray and 7- 
ray light curves can be reproduced, given the inclination 
angle is not large. In contrast, the OG model (Cheng et 
al. 1986a, CHR) presumed that the gaps can exist in the 
outer magnetosphere between the null charge surface and 
the light cylinder. Later, CHR model was developed to the 
single OG models (e.g., Chiang & Romani 1994, Romani & 
Yadigaroglu 1995, Cheng et al. 2000), which claimed that 
the three dimension extents of OGs are constrained by the 
pair cascade processes, and a single OG can produce wide 
double-peak high energy light curves. Whereas theoretical 
considerations of more detailed physical processes for par- 
ticle acceleration and photon emission are necessary, it is 
urgent and interesting to find new observational evidence 
for these models. 

PSR B1937-I-21, with the period of 1.56 ms, is the fastest 
millisecond pulsar (MSP) known. At radio bands, it ex- 
hibits an interpulse emission that is roughly equal to the 



main pulse in intensity and separated from it by a phase 
of about 180°. Recently, the ASCA observations detected 
pulsed X-ray emission from this pulsar (Takahashi et al. 
2001). Only one nonthermal narrow pulse was observed, 
which is coincident with the radio interpulse in phase 
within the timing errors. The pulse width is about 100±61 
(23° ± 14°). Besides the narrow pulse, the light curve 
reveals two additional wide Gaussian-shaped bulges above 
the background level, with each phase interval being about 
0.5 rotation period and each peak intensity ~ 1/4 of the 
nonthermal pulse peak. 

Where does the X-ray emission of the fastest MSP come 
from, the PC or the OG? It has not been extensively stud- 
ied as the other high energy pulsars such as the Crab pulsar 
and Vela. The previous work was done by Luo et al (2000), 
in which the PC model was modified for MSPs. They ap- 
plied the theory to PSR B1937-f 21 and suggested that the 
X-ray emission probably originates from the location of 
one stellar radius above the PC. 

In this paper, efforts on modelling the observational 
data at radio and X-ray bands are made for both the PC 
and OG models. The inclination angle between the rota- 
tion and magnetic axes is constrained in section 2, which 
is necessary for the modelling. The modelling is carried 
out in section 3. Conclusions and discussions are placed 
in section 4. 

2. THE INCLINATION ANGLE 

The inclination angle (a) between the rotation and mag- 
netic axes is a necessary parameter for both the PC and 
OG models to give various high energy emission beams. 
Unfortunately, there is no agreement on the value of incli- 
nation angle of PSR B1937-I-21. In this section we rein- 
vestigate the value of a under the double-pole model, viz., 
the radio interpulse and main pulse are considered to be 
from the opposite magnetic poles of a dipolar field, and 
the result is used in the calculation in section 3. 

We assume that the radio emission beams from double 
poles are axisymmetric around the magnetic axis and have 
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the same radius, namely, pi = p2- According to the ge- 
ometry model (Gil et al. 1984, Lyne & Manchester 1988, 
LM88) one has 
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where /3 is impact angle between the line of sight (LOS) 
and the magnetic axis, W is the pulse width of the aver- 
age profile, the subscripts '1' and '2' denote the main pulse 
and the interpulse, respectively. There are two simple ge- 
ometrical relations between the inclination angles and the 
impact angles, namely, a2 = tt — ai, and /32 = Q!i-|-/3i — 02- 
^Prom the above relations (3i can be derived as 
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In the following we neglect the subscript '1', so that all 
the a, (5 and p below are referred to the main pulse except 
for special declaration. 

To figure out the value of A, the profile at 1.5 GHz is 
chosen (downloaded from EPN) for its high time resolution 
and low dispersion smearing (Kramer et al. 1998). The 
pulse widths arc mcasrued at the level of 10% of their peak 
intensities, which are 19°.4 ± 0°.4 and 17°.8 ± 0°.4 for the 
main pulse and the interpulse respectively. So that /3(a) 
and p{a) can be calculated, as shown by the solid curve 
and the curve dotted by circles respectively in Fig.l. The 
(3{a) curve approximates to a linear function of a+f3 « 90° 
when a < 85°, which is determined by the fact that the 
main pulse and interpulse are both narrow and differ only 
a little in pulse width. 

For a dipolar field, the shape of the polar cap is found 
to change from a circle (for a = 0°) to an ellipse of which 
the longitudinal radius is about 1.6 times of the latitudi- 
nal radius (for a = 90°, Cheng et al. 2000). Since the 
deviation from circular shape is not essentially significant, 
we simply regard the polar cap as a circle in this section. 
So that the opening angle (between the magnetic axis and 
the tangent of the magnetic field line) of the polar cap ppc 
can be determined by the last open field line on the plane 
containing the rotation and magnetic axes {Q, — /i plane) . 
The radius ppc is a function of the stellar radius R and 
the inclination angle a, as shown by the dash curves in 
Fig.l for i? = 3 km and 10 km. 

We further assume that the boundary of the radio beam 
is defined by the last open field lines, then the beam radius 
should be greater than ppc- According to this criteria, the 
inclination angle is constrained to be a < 63° for i? = 10 
km and a < 77° for i? = 3 km. As the radii of neutron 
stars (NSs) are currently believed to be about 10 km, we 
accept that a < 63°. 

Can we determine the exact value of al It is sure if 
the maximum rate of position angle (dV'/d(/))max is exactly 
known, where ^ is the position angle and (j) the pulse longi- 
tude. The maximum rate presents the second relationship 
between a and /?, which reads 

'di/^X sin a 



Combining with Eq.s 3 and 5, a and (3 can be solved for 
a given (di/'/d0),„ax- 

Although recent polarization observations present flat 
position angle sweeps (Thorsett & Stinebring 1990, Stairs 
et al. 1999), it can not be asserted that the real value 
of (dV'/d^)max is small, because observations may give a 
less steep position angle gradient due to smearing of finite 
sampling time, to the frequency dispersion in pulse arrival 
time (Liu & Wu, 1999), and to the interstellar scatter- 
ing (Gil 1985a). The /3'(q:) cm-vcs derived from Eq. 5 for 
(d'!/'/d(/>)inax = 1,3 and 20 respectively are presented by 
the dotted curves in Fig.l. The intersections of /3'(a) and 
j3{a) shows that a larger (d'(/'/d(/))max results in a larger 
a and smaller [3, which means the LOS sweeps across the 
radio beam more closely to the bc;ani ccinter. 

One may find that when the real value of (dV'/d^)max 
is large enough, a would exceed 63° (for example, taken 
(d'(/'/d(/))max = 3 as proposed by Gil (1985a), a is 71°) 
and hence contradict against a < 63°. However, if the ra- 
dius is smaller, for example, i? = 3 km, this inconsistency 
would cancel. In fact, Xu et al. (2001) suggested that 
PSR B1937-F21 is probably a strange star (SS) with low 
mass and small radius. The detailed discussion is placed 
in section 4. 

The range of a presented above is different from the con- 
ventional consideration in double-pole model that a should 
be close to 90° (Stairs et al. 1999). Alternatively, there 
is anotlica' kind of so-called single pole; model to interpret 
the interpulse, which suggests that the interpulse emission 
comes from the same pole as the main pulse. In the single 
pole model proposed by Gil (1985a) for PSR B1937+21, 
a only need to be 20°. The single pole model predicts 
that the separation between the main pulse and interpulse 
may be close to 180° and is frequency independent (Gil 
1983, 1985b), which is coincident with the observation of 
PSR B1937+21 (Hankins & Fowler 1986). However, ob- 
servations with high time resolution (e.g. Kramer et al. 
1998, Stairs et al. 1999) failed to find the emission com- 
ponents between the main pulse and interpulse, which was 
reported by Stinebring et al. (1984) and was suggested to 
be a strong support to the single pole model (Gil 1985a). 
Therefore, in this paper the radio emission is considered 
to come from double poles. More confirmative estimations 
of (dV'/d(^)inax are expected to determine a and (3. 



3. THE ORIGIN OF NONTHERMAL X-RAY EMISSION FROM 

PSR B1937+21 

In this section we calculate the X-ray beams in the 
frames of both the PC and OG models to find out whether 
they are able to reproduce the observational facts. The 
facts used here are: 

(1) at 1.4 GHz the separation between the peaks of in- 
terpulse and main pulse is 174° (measured from the profile 
presented by Takahashi et al. (2001)); 

(2) at 1.5 GHz (EPN data) the 10% widths of the main 
pulse and interpulse are 19.4° and 17.7°, respectively; 

(3) the nonthermal X-ray pulse is nearly coincident with 
the interpulse, the X-ray pulse width is about 23° (Taka- 
hashi et al. 2001). 
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3.1. Origin from the PC accelerators? 

Luo et al. (2000) discussed the viability of PC models 
for high energy emission from MSPs. They found that the 
maximum Lorenz factor of particles is limited by curva- 
ture radiation and not sensitive to the specific acceleration 
model, but the height where the Lorenz factor achieves 
the maximum is model dependent, which may be between 
O.Oli? (for the inner vacuum gap) and above O.li? (for the 
space-charge limited gap) from PC surface for pulsar pe- 
riod P = 2 ms and a surface magnetic field Bg = 7.5 x 10^ 
G. Assuming a space-charge limited flow, the pair cas- 
cades can occur at the typical distance (to the star center) 
of r ~ (1.5 — 2.5)i?, and high energy emission is radi- 
ated from this region. With respect to P = 1.56 ms and 
= 4.1 X 10^ G for PSR B1937-h21, their analysis applies 
to this pulsar. 

Since the radio emission is also radiated from the re- 
gion near the PC, homocentric radio and X-ray beams are 
produced (Fig. 2a). Quantitatively, we assume the radio 
emission arises from the PC up to the distance r = 2.5P, 
the X-rays may arise from the PC to a (a) relatively higher 
distance, e.g., r = 3R, or (b) lower distance, e.g., r = 2P. 
For simplicity we assume the emission regions are bounded 
by the last open field lines. In case (a) the X-ray beams 
are wider than the radio beams, therefore, when the line 
of sight sweeps across the both radio beams (to reproduce 
the main pulse and the interpulse) it would sweep across 
the both X-ray beams either, and hence gives double X-ray 
pulse, which is inconsistent with the observation. In case 
(b) the X-ray beams are narrower than the radio beams, 
and then, for a proper viewing angle (between LOS and 
the rotation axis), only one X-ray beam could be observed. 
In this case, could the fact (3) be explained as well? We 
analyze this issue by using the figure of the (phase- viewing 
angle) plane on which the emission beams are projected. 

In Fig. 3, the horizon axis is pulse phase, the vertical 
axis is is the viewing angle ( {( = a + f3). The dotted 
curves represent the boundaries of X-ray beams, the solid 
ciirvcs are the boundaries of radio beams, and the PCs are 
shown by the dashed curves. Retardation (due to distinct 
emission heights) and abberation effects (due to emission 
sources co-rotating with the pulsar) have been taken into 
account, both of which make the lower beam move to- 
wards the trailing edge of the higher beam. To obtain the 
figure, a moderate inclination angle is taken, a = 60°. Al- 
though there is a < 63° as discussed in section 2, the angle 
should not be too small, or it would give too wide radio 
beams (Fig.l), which must be emitted from unbelievably 
high distances near the light cylinder. 

According to Fig. 3, when ( ~ 80° the X-ray peak is 
coincident with the radio pulse centered on phase about 
0.9; when C ~ 100° it is coincident with the radio pulse on 
phase about 0.4. However in both cases, the radio pulse 
associated with the X-ray pulse is the main pulse which is 
wider and more intense than the other one. 

However, we should point out that the above analysis 
is based on a simple assumption that the radio (X-ray) 
beams from both poles have the same width and the emis- 
sion pattern in the beams is symmetrical about the mag- 
netic axis. Since the magnetosphere structure of pulsars 
and the detailed radio emission process are still unknown 
exactly, possibly significant deviation from the assumption 



can not be ruled out, and the validity of the PC model 
need to be considered further. We discuss such possibility 
in section 4. 

3.2. Origin from the OG accelerators 

In the original OG model CHR suggested that a global 
current flow through the magnetosphere can result in large 
regions of OGs between the null charge surface and the 
light cylinder along the last open field lines. Within the 
OGs particles with one kind of charge are accelerated out- 
ward from the star and give an outward emission beam, 
while those with the opposite charge are accelerated to- 
wards the star and give an inward beam. The high energy 
photons were proposed to be emitted from two OGs associ- 
ated with the two poles so that double-peaked 7-ray pulse 
profile can be reproduced, of which one peak corresponds 
to the outward beam from one OG and the other peak 
to the inward beam from the opposite OG. CHR assumed 
that the OGs are active only near the ft — jj. plane. How- 
ever, this assumption is merely valid for large inclination 
angles. 

The currently prevalent OG models are the single OG 
models (e.g. Chiang & Romani 1992, Chiang & Romani 
1994, Cheng et al. 2000). Generally the inward emission 
is not important in these models for the reason that the 
inward high energy photons can not pass freely through 
the inner magnetosphere due to magnetic pair production. 
The outward emission from the OG associated with a sin- 
gle pole can produce a broad, irregularly-shaped emission 
beam which is particularly dense near the edge. The OG 
regions can be supported along all the last open field lines, 
but the three dimension scales of OGs are limited by the 
pair production processes. 

In the latest version of this type of model (Cheng et 
al. 2000), the fraction size (/ s h/Ri,) of the gap is 
/ ~ 5.5P26/2iBj"//^C"i/^ which can be estimated by the 

threshold of 7— 7 pair production, i?x(/)-£'7(/) 5^ (n^cC^)'^, 
where h is the mean vertical extension perpendicular to 
the magnetic field, Pl is the radius of light cylinder, 
^ = A(/)/27r, A(t> is the transverse extension of the gap, 
is the energy of the X-ray photons emitted from hot PCs, 
and Ej is the characteristic photon energy emitted by the 
relativistic particles. The radial scale of pair production 
is limited to a range between rin and rum ~ 6rin(^ = 0), 
where rjn (the siibscript 'in' means the inner boundary 
of the OG) is the distance of null charge surface, ^ = 
corresponds to the ft — fj, plane. 

In the following our modelling is in the frame of single 
OG model. Only the outward emission beams from two 
OGs are considered, as illustrated by Fig. 2b. For PSR 
B1937+21, we have a thin OG with'/ = 0.16 ^^^Z^, so 
that the X-rays can be simply regarded as being radiated 
from the last open field lines unless ^ is too small. The 
radial scale rii,„/rin and the transverse scale ^ are free pa- 
rameters in calculating the X-ray beams. 

First we consider a general situation of the OG scenario 
to test if the observation facts can be hopefully repro- 
duced. From Fig. 2b one can see that provided the ob- 
server's viewing angle is not just 90°, e.g., C = 83°, the 
LOS can sweep across both of the radio beams and only 
one X-ray beam. By assuming a group of reasonable val- 
ues of parameters (see the first line of Table 1), namely, 
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the inclination angle, the stellar radius, the distance of 
radio emission, and the radial and transverse extensions 
of the OGs, the X-ray beams are calculated, which are 
demonstrated by the line-shadowed areas in Fig.3. Retar- 
dation and abbcration effects are also included. It shows 
clearly that the X-ray pulse could be associated with the 
radio interpulse (which has a smaller width than the other 
one), and the narrow X-ray pulse width may be obtained 
if proper extension of the OG is assumed. 

Then wc model the observational data. The parameters 
listed in the second line of Table 1 are found out to be able 
to reproduce the narrow X-ray pulse width, the radio pulse 
widths, and the coincidence between the X-ray pulse and 
the radio interpulse, which arc in good agreement with the 
observational data. In order to simulate the observational 
profiles, we simply assume that the X-ray and radio pulses 
are Gaussian shapes, and assume additionally wide, weak, 
hot X-ray emission from both of the PCs. The theoretical 
profiles are plotted in Fig. 4, together with the observa- 
tional profiles for comparison. 

It should be pointed out that modelling the X-ray pulse 
width is not sensitive to the value of riim/rjn but to ^, 
thereby the range of riim/rin is relaxed and a reason- 
able value is chosen. Other groups of values are also 
tried. It is found that for the moderate inclination angles 
40° < a < 63° the observational data can be reproduced 
with proper choice of the gap size = 40° in Table 1 
is approximately the maximum transverse scale). There- 
fore, our calculation suggests that the nonthermal X-rays 
of PSR B1937-F21 may be emitted from the OGs. 

4. CONCLUSIONS AND DISCUSSIONS 

The discovery of pulsed X-ray emission from the MSP 
PSR B1937-I-21 and the phase alignment between the X- 
ray pulse and the radio interpulse provide valuable infor- 
mation on the emission mechanism. In this paper we in- 
vestigate the possible origin of the X-rays from both the 
PC and OG accelerators. In the frame of the prevailing 
OG model (e.g. Cheng et al. 2000), the X-rays from outer 
gap accelerators could account for the main observational 
facts by assuming proper size of the OGs for an oblique 
rotator. In the frame of PC model, by assuming sym- 
metric geometry for radio and X-ray emission, the X-rays 
from polar cap accelerators would produce an X-ray peak 
aligned with the main radio pulse in phase, which is con- 
tradictive against the observation. 

Some more discussions for both the OG and PC models 
are presented as follows. First we refer to the OG model. 

There is a slight inconsistency as shown in Fig. 4, i.e., 
in the OG model, the calculated separation between the 
interpulse and main pulse is 180°, 6° greater than the ob- 
servational value. This may be due to the retardation 
effect. The phase shift for different heights can be roughly 
estimated as As = Ar/ (Pc). A difference of Ar = 0.8R is 
enough to producx; the phase shift of 6°. 

A moderate inclination angle 55° is used in the mod- 
elling. But a could be larger if the real gradient of position 
angle is steeper than the present observations as suggested 
in section 2. Could the OG model be still valid for large 
value of a? 

As shown in Fig.l, when (d'!/'/d(/>)inax ^ 4 it would result 
in a puzzling problem that the derived radio beam radius 



is considerably smaller than the PC radius, if i? = 10 km 
and a magnetic dipole are assumed. A much smaller stellar 
radius could cancel the problem, but this requires an SS 
scenario, because the smallest radius of NSs is generally 
believed to be ^9 km while SSs can have much smaller 
radius due to their different equations of state. In fact, 
according to the observational limits on the radius and 
mass derived from the pulse width and (dV'/d^)max, Xu 
et al. (2001) suggested that PSR B1937-h21 is probably a 
strange star with much low mass, small radius and weak 
magnetic moment. If PSR B1937-I-21 is an SS with small 
radius, for example, R = 2 km, a group of parameters is 
found out to be able to model the observational data by 
OG model, which is listed in the third line in Table 1. 

Then we turn to the PC model. When the emission ge- 
ometry is not symmetrical, the PC model could be able to 
explain the observation as well. The possible asymmetry 
is discussed for the radio and the X-ray emission, respec- 
tively. 

(1) It has been assumed that the radio emission pat- 
tern from the two poles is the same in our modelling; the 
different behaviors of the main pulse and the interpulse 
are thus geometrical origin. However, their difference may 
be caused intrinsically since the mechanism for the radio 
emission is still poorly understood. For example, a possi- 
ble reason may be that the pulsar has a non-dipolar mag- 
netic field, thus the distribution of emission intensity in 
the radio beams may be different from each other, which 
could lead to that for one beam only part of it is observed 
while for the other a larger part or the whole is observed. 

(2) When PSR B1937-h21 is an NS or an SS with crust 
and its inclination angle is near 90°, the X-ray emission 
pattern may be significantly asymmetric about the mag- 
netic axis. In this case, the space charge is negative on the 
side toward the rotation axis (where n • B > 0, hereafter 
side I) and positive on the side away from the axis (where 
r2 • B < 0, hereafter side II). On side II the ions can be 
pulled away from the surface by strong electric field if the 
binding energy is small enough. Therefore, the half beam 
from side II could be much less luminous than that from 
side I due to the much smaller Lorentz factors of the ions, 
and thus may be too weak to be observed. Notice that, if 
side I is above the equator on one pole, it should be be- 
low the equator on the other pole. Therefore, only single 
X-ray peak is observed which may coincide with the radio 
interpulse. 

In recent years, it is suggested that pulsars may be bare 
strange stars (BSSs) (e.g., Xu 2002 and references therein). 
If PSR B1937+21 is also a BSS, the positive charge on side 
II is carried by positrons but not ions, then the emission 
pattern on both sides should be the same, and each beam 
may be symmetric around the magnetic axis. In this case, 
when symmetrical geometry is assumed for the radio emis- 
sion, PC model could not account for the observational 
facts, otherwise, PC model may still be valid. 

In general, further research on the pulsar's magneto- 
sphere structure and emission mechanisms will be helpful 
to understand the origin of its X-rays for PSR B1937-I-21. 
We are expecting that future polarization observations 
could provide more confirmative value of {dip/d(p)aiax, 
which is meaningful not only for constraining whether this 
pulsar is an NS or an SS but also for a better understand- 
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ing of its X-ray emission. 
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Table 1 

Parameters for calculating the radio and X-ray beams/profiles 



No. 


an 


C(°) 


R (km) 


(i?) 


nn(0)'^(i?) 




(°) 


1 


60 




10 


2.5 


1.3 


3.5 


100 


2 


55 


89.5 


10 


1.7 


1.5 


2.0 


40 


3 


75 


89.5 


2 


2.9 


2.1 


2.5 


66 



''the distance of radio emission, in unit of stellar radius R. 
^the distance of the null charge surface on the 17 — /x plane. 
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Fig. 2. — (a) Scheme for the X-ray beams in the frame of PC model, (b) Scheme for the X-ray beams produced by the OG model. Radio 
beams are also plotted. 
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Fig. 3. — Emission beams projected onto the (phase-viewing angle) plane for a = 60°. The solid curves are the boundaries of radio beams, 
the dotted curves are the boundaries of X-ray beams from the extended acceleration zones above the PCs suggested by the PC model, the 
line-shadowed areas are the X-ray beams from the outer gaps, and the dashed curves represent the PCs (See text and the first line of Table 
1 for the related parameters). 
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Fig. 4. — Theoretical X-ray and radio profiles (the dashed curves) together with the observational profiles (Takahashi et al. 2001). The 
parameters in the second line of Table 1 are used to obtain the theoretical profiles. 



